Geometrical derivation of the Boltzmann factor 
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Abstract 

We show that the Boltzmann factor has a geometrical origin. Its derivation follows from the 
microcanonical picture. The Maxwell-Boltzmann distribution or the wealth distribution in human 
society are some direct applications of this new interpretation. 
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The Maxwellian distribution was derived in Ref. Q by geometrical arguments over an A^- 
sphere. Following this insight, here we explain the origin of the Boltzmann factor by means 
of the geometrical properties of an A^-hyper plane. This alternative view can be thought as an 
a priori theoretical argument that explains why the probability of a magnitude, constrained 
by a conservation law, is given by the exponential distribution when the ensemble or system 
over which that magnitude is measured is in thermodinamical equilibrium. 

By thermodinamical equilibrium we mean the macroscopic stationary regime where all 
the microscopic accessible states of the system are equiprobable. This implies that the 
deterministic or random time evolution of the system verifies the ergodic hypothesis. This 
has been proved in very few systems^, and hence, in general, it is supposed as an hypothesis 
in classical statistical mechanics-. If the system do not explore its own phase space with 
equiprobable results for all the accessible states then the system exhibits some kind of 
ergodicity breaking and the final distribution will not be the exponential distribution. We 
can call it a non- Boltzmann equilibrium. 

Without loss of generality, here we obtain analytically the Boltzmann factor (or Maxwell- 
Boltzmann distribution) in a model where the ergodic hypothesis has been computationally 
proved. The method was developed in Ref. !li and it is grounded on geometrical arguments. 

The model was proposed by Dragulescu and Yakovenko^ in order to explain the distribu- 
tion of wealth in human society. Thus it is known that the incomes of 90% of the population 
in western societies can be fitted by a exponential distribution.- Supposing initially equity in 
a set of many agents, that is, all of them have the same initial quantity of money, the model 
establishes random binary interactions in which agents exchange money but conserve its 
total amount. It is computationally found that the system asymptotically tends toward the 
exponential distribution. This means that, in this case, the evolution mechanisms proposed 
by Dragulescu and Yakovenko for the money exchange among agents verify the ergodic hy- 
pothesis and, then, its final distribution is the exponential distribution. (Let us say that 
other non-Boltzmann random mechanisms for the money exchange are also proposed in 
Ref. 4- These mechanisms imply the breaking of the ergodic hypothesis and therefore the 
asymptotic evolution of the system toward non-Boltzmann equilibria). 

We start by assuming agents, each one with coordinate (money) Xi, i = 1, . . . ,N, with 
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Xi > 0, and a total amount of money E that is conserved: 



Xi+X2-\ \-Xn-1+Xn = E. (1) 



Under the evolution rules proposed in Ref. |4, this isolated system evolves on the positive 
part of an equilateral A^-hyperplane (i.e., the surface formed for all those points with Xi > 0, 
for all i). The formula for the surface area Si\i{E) of an equilateral A^-hyperplane of side E 
is 



(See appendix A for the derivation of this formula). If the ergodic hypothesis is assumed, 
each state of the microcanonical ensemble (that is, each point on the A^-hyperplane) is 
equiprobable, then the probability f{xi)dxi of finding the agent i with money Xi is pro- 
portional to the surface area formed by all the points on the A^-hyperplane having the 
zth-coordinate equal to Xj. Our objective is to show that /(xj) is the Boltzmann factor (or 
the Maxwell-Bolztamnn distribution), with the normalization condition 

rE 



f{xi)dxi = 1. (3) 

If the ith agent has coordinate Xi, the — 1 remaining agents share the money E — Xi 
on the {N — l)-hyperplane 

X1+X2 h Xi-i + Xi+i ^ xn = E - Xi, (4) 

whose surface area is Sn-i{E — Xi). If we define the coordinate 6]^ (see appendix A) as 
satisfying 

In -I 

sin^jv = y^^, (5) 

it can be easily proved that 

S^{E) = rSj,^,{E-x,)^. (6) 
JO smfc'iv 

Hence, the surface area of the A^-hyperplane for which the ith coordinate is between Xi 
and Xi + dxi is Sn~i{E — Xi)dxi/ sin 6'Ar. We rewrite the surface area as a function of Xj, 
normalize it to satisfy Eq. ([3]), and obtain 

r, ^ _ 1 SN-l{,E-Xi) 



whose final form, after some calculation is 



fix.) = iN-l)E-'{l-^f-\ (8) 
If we call e the mean wealth per agent, E = Ne, then in the limit of large we have 

/ ™.\Af-2 

limfl--) ~e-"^/^ (9) 

The Boltzmann factor e~^*/^ is found when ^ 1 but, even for small A^, it can be a good 
approximation for agents with low wealth. After substituting Eq. into Eq. ([8]), we obtain 
the Maxwell-Boltzmann distribution in the asymptotic regime N ^ oo (which also implies 
E oo): 

f{x)dx=-^e~''''dx, (10) 

where the index i has been removed because the distribution is the same for each agent, and 
thus the wealth distribution can be obtained by averaging over all the agents. 

Depending on the physical situation the mean wealth per agent e takes different ex- 
pressions and interpretations. For instance, in the case of an isolated one-dimensional 
gas, the Xi variable represents the energy of the particle z, thus we can calculate the de- 
pendence of e on the temperature, which in the microcanonical ensemble is defined by 
differentiating the entropy with respect to the energy. The entropy can be written as 
S = — kNJ^^ f{x) In f{x) dx, where f{x) is given by Eq. fIlOp and k is the Boltzmann 
constant. If we recall that e = E/N, we obtain 

S{E) = kN In (^)+kN. (11) 
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The calculation of the temperature T gives 



Thus e = kT, and the Maxwell-Boltzmann distribution is obtained as it is usually given in 
the literature: 

f{x)dx=^e-^/'''dx. (13) 

In conclusion, we have shown that the Boltzmann factor describes the general statistical 
behavior of each small part of a multi-component system whose components or parts are 
given by a set of random variables that satisfy a conservation law (condition ([1])), and that 
run in an equiprobable manner over all the states of its own phase space. 



APPENDIX A: Formula for the surface of an equilateral A^-hyperplane 



Here we derive the formula for the calculation of the surface Sn{E) of an equilateral A^- 
hyperplane, Ujy, embedded in and formed by all the points {xi,X2, ■ ■ ■ ,X]^), with Xj > 
for alH = 1, 2, . . . , A^, verifying the equation 

Xi + X2-\ \-Xn-1+Xn = E, (1) 

with E a real constant. 

First we define the angle dj^i as follows. The unitary perpendicular vector to the hy- 
perplane Un in 3?^, u±, has coordinates 1, • • • , 1, !)• The unitary vector in the xn- 

direction, Vxj^, has coordinates (0, 0, ... , 0, 1). The projection of lj± over t;^^ is given by the 
scalar product of both vectors: 

cos On = uj±- v^^ = , (2) 



that defines the angle 9n- Let us observe the curious property of this angle: 

cos^^ . . 
= sm6'Ar. (3) 

COSt^Ar_l 

The surface Si\f{E) in is the measure of a set of points with — 1 dimensions, just 
the set of all those points {xi,X2, ■ ■ ■ ,xn), with Xi > 0, verifying Eq. ([T]). The one-to- 
one projection of this set of points on the XAr-direction generates a set of points of — 1 
dimensions in 3?^"^ just the set of all those points in 3?^~^ satisfying the inequality: 

X1 + X2 + \-xn-i<E. (4) 

The volume of this set of points Vn~i{E) is the volume of an (A^ — l)-dimensional pyramid 
formed by A^ vertices linked by A^ — 1 perpendicular sides of length E. It is straightforward 
to see that the formula for this volume is: 

T?N-1 

Moreover, Vn^i{E) and Sn{E) can be calculated explicitly when A^ = 1, 2, 3, 4. They verify 
the relationship: 

VN-i{E) = SN{E)-coseN. (6) 
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This expresses an exact connection between the surface Sn{E) of the A^-hyperplane in 
and its projection in the XAr-direction over the (A^ — l)-dimensional volume Vn^i{E) in 3?^^^ 
for those low dimensional cases. Extrapolating this behavior for all N, we conclude that 

as it has been used in Eq. ([2]) of the main text. 

Although we can not claim the mathematical proof of formula ([7]), the Boltzmann factor 
is exactly obtained in the limit ^ 1 after taking this formula as starting point of all 
the line of reasoning developed in the main text. This result can be used like a posteriori 
argument that proves the correctness of expression (I7j). 



APPENDIX B: A possible generalization: an open problem 

In this work, we have shown that an ensemble of positive variables (xi, X2, . . . , xn) veri- 
fying 

Xl + X2^ \-Xn-1+Xn = E, (1) 

with an adequate mechanism assuring the ergodic hypothesis, i.e., the equiprobability of 
all the possible states (xi,a:2, . . . ,xn) on the hypersurface of the phase space, presents the 
exponential distribution of the generic variable x over the ensemble when N oo, 

f{x)dx ~ e"^ e-^/^rfx, (2) 

with E = Ne and e the mean value of the real constant E over the collectivity. 

In the work with Ref. Q, it was shown that an ensemble of positive variables 
{xi,X2, . . .,xn) verifying 

Xi -\- X2 ~l~ ■ ■ ■ ~l~ Xj^_i -\- X = E^ (3) 

with an adequate mechanism assuring the ergodic hypothesis, i.e., the equiprobability of 
all the possible states (xi,X2, . . . ,xn) on the hypersurface of the phase space, presents the 
Gaussian distribution of the generic variable x over the ensemble when N 00, 

f{x)dx ~ e-y^^dx, (4) 

with E = Ne and e the mean value of the real constant E over the collectivity. 



6 



The general question that we want to bring to the reader is the following. Let 6 be a real 
constant. If we have a set of positive variables (xi, X2, ■ ■ ■ , xn) verifying 

x^ -\- X2 ~l~ ■ ■ ■ ~l~ Xj\f_-^ -\- X j\j = E (5) 

with an adequate mechanism assuring the ergodic hypothesis, i.e., the equiprobability of all 
the possible states (xi,X2, . . . ,X]^) on the hypersurface of the phase space, will we have for 
the generic variable x the distribution 

/(x)dx~e-^/'e-"'/'^da;, (6) 

when we average over the ensemble in the limit 00?. The answer to this question, as 

far as we know, is an open problem. 
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